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We present an extension of the c-axis tight-binding sum rule discussed by Chakravarty, Kee, 
and Abrahams [Phys. Rev. Lett. 82, 2366 (1999)] that apphes to multilayer high-Tc cuprate su- 
perconductors (HTCS) and use it to estimate — from available infrared data — the change below Tc 
of the c-axis kinetic energy, (He), in YBa2Cu307_4 (5 = 0.45,0.25,0.07), Bi2Sr2CaCu2 0g, and 
Bi2Sr2Ca2Cu30io. In all these compounds (He) decreases below Tc and except for Bi2Sr2CaCu208 
the change of (He) is of the same order of magnitude as the condensation energy. This observation 
supports the hypothesis that in multilayer HTCS superconductivity is considerably amplified by the 
interlayer tunnelling mechanism. 



There is a growing evidence that at least in some HTCS the effective out-of-plane (c-axis) kinetic energy, (He), 
decreases upon entering the superconducting (SC) state In this manuscript we address the question whether 

this decrease represents a significant (and perhaps the dominant) contribution to the condensation energy or whether 
it is merely a small byproduct of an in-plane pairing mechanism. It is obvious that a quantitative estimate of the 
change of (He) between the normal and the SC state, A{Hc) = {Hc){T w Tc) — {Hc){T « Tc), is required in order 
to answer this interesting question. 

Two different approaches have been recently employed to obtain the value of A{Hc)'. an approximate one Q where 
A{Hc) is identified with the Josephson coupling energy ("JCE", Ej) of the internal Josephson junctions and a rigorous 
one based on the so called c-axis tight-binding sum rule [Q. Within the JCE approach, the value of A{Hc) per unit 
cell is given by the formula 

MHc) = E.J = "^ul (1) 

which contains only one nontrivial parameter, the plasma frequency of the internal Josephson plasmon ojpi. The values 
of the in-plane lattice constant a and the distance between the neighboring copper-oxygen planes d are well known. 
The sum-rule (SR) approach instead relates A{Hc) to the increase below Tc of the low- frequency optical spectral 
weight (SW) 

nijj 

a{T,cu) = / aic{T,iu')dio' , (2) 



Jo 

A{Hc) = ^^HT « Tc,nc)-a{T « Tcilc)] = ^^Aa. (3) 

, Here flc is a cutoff frequency and aic is the real part of the c-axis conductivity cTc- What is the connection between 
^ the two approaches? The JCE approach yields only the contribution of the condensate (related to the spectral weight 
5^ Ps of the 6 peak at w = in aic^Lo)) but neglects the contribution due to the single particle tunnelling that is related 
to the change AiV of the finite frequency part N{T, fie) of a{T, fie)- The SR approach, on the other hand, considers 
both contributions to A{Hc). Only in the limit of vanishingly small changes at Tc of the regular part of Cc the two 
approaches can be expected to yield the same result (except for a factor of 4 as pointed out in Ref. Q). Especially 
the underdoped cuprate HTCS are not very far from this limit where the JCE approach can be expected to yield a 
reasonable — better than order of magnitude — estimate of AHc- 

In a previous paper |^ we and our coworkers have reported the values of Ej of two compounds that have two 
copper-oxygen planes per unit cell (bilayer compounds): YBa2Cu307-5 (Y-123) and Bi2Sr2CaCu208 (Bi-2212). Note 
that in these compounds Ej is determined mainly by the frequency uiu of the intra-bilayer Josephson plasmon |q-pT|] 
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and by the distance between the closely spaced copper-oxygen planes. We have shown that there is a reasonably 
good agreement between the values of Ej and the values of the condensation energy Uq obtained from the specific 
heat data. Assuming A{Hc) ~ Ej we arrive at the conclusion that the condensation energy (and the high value of 
Tc) in the two compounds can be accounted for by the change at Tc of {He), i.e., by the interlayer tunneling theory 
| fi2[ . Nevertheless, it can be objected that the actual value of A{Hc) may be much smaller than the one of Ej and 
therefore should be determined by using the complementary and more rigorous SR approach. Unfortunately the 
approach in its simplest form as presented above does not apply to multilayer compounds |^ . Its derivation is based 
on the assumption that the distribution of the total internal electric field in the superconductor is homogeneous, a 
condition that is not met for multilayer systems. In the following we present an extension of the SR approach that 
allows one to obtain a reasonable estimate of A{Hc) in multilayer cuprate HTCS. We provide a short derivation of 
the key formulas and discuss the applications to Y-123, Bi-2212, and Bi2Sr2Ca2Cu30io (Bi-2223). 

For the sake of simplicity we focus on the bilayer case and we use the notation introduced in Rcf . Q , in particular 
in its appendix B. The effective c-axis kinetic energy (He) of a bilayer compound consists of the intra-bilayer term 
{Hbi) and the inter-bilayer one {Hint), 

H, = Hm + H,nt and A{H,) = A{Hm) + A{H„u) ■ (4) 

The current densities in the intra-bilayer region and in the inter-bilayer region denoted as jbi and jint , respectively, are 
related to the two electric fields Eu and Eint as follows (Eq. (B7) of Ref. |^): ja{i^) = (Tap{uj)Ef3{cv), a, /? e {bl, int}. 
It is reasonable to assume that the conductivity matrix ctq (3 is diagonal and to introduce the abbreviations au = cFbi u 
( "intra-bilayer conductivity" ) and (Tint = Ointint ( "inter-bilayer conductivity" ) so that jf,; = auEu and ji„t = UintEint- 
The total conductivity is then given as follows: 



a{uj) = {du + dint) I 



dbl ^ dint 
(Tbl{^) 0-„t(w) 



(5) 



which is the formula used at the phenomenological level introduced in Ref. ||^ . We aim at finding a formula connecting 
A{Hc) and Aa as defined in Eq. (3). Let us first express the two components A{Hbi) and A{Hint) of A{Hc) in terms 
of the two conductivities. It follows from Eqs. of Appendix B of Ref. Q that 

A{Hbi) = ^^^Aau (6a) 

and 

A{H,nt) = AaM , (6b) 

where cxbi/int and Aabi/int are related to cri,i/int{^) iii the same way as a and Aa are to ac{uj). A schematic 
representation of Aa;,; and Aaint is shown in Fig. 1. The quantity Aa can also be expressed in terms of Aa;,; and 
Aaint- After some manipulations (similar to those in chapter 5.7 of Ref. [l3| ) using the analytic and the asymptotic 
properties of cr(w), abi^to) and CTmt(w) we obtain 

a{T,uj —^oo)= — abi{T,uj ^ co) + ^"""^ — ai„i(r, w oo) . 

dbl + dint dbl + dint 

Provided that the temperature dependence of the three conductivities above fic is negligible we can also write 

= ■ ^''j — + , — Aaint ■ (7) 

dbl + dint dbl + dint 

In order to obtain an estimate of A{Hc) we have to make an additional assumption because the Eqs. (4), (6a), (6b) 
and (7) still contain five unknowns: A{Hc), A{Hbi), A{Hint), Aabi and Aaint {Aa is assumed to be known from the 
infrared data). We suggest to use the following assumption: 



Aaint 



(8) 



int 



where lou and tOint are the superfluid plasma frequencies of the intra-bilayer region and the inter-bilayer region, 
respectively. In other words, we suggest that Aa;,; /int is proportional to the contribution i,i jint of the condensate 
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(see Fig. 1 for a definition of Psu/int)- Using Eqs. (4), (6a), (6b), (7) and (8) we finally arrive at the formula expressing 
A{Hc) in terms of Aa: 

dbldint ^bldbl + ^intdint 

where k = 27i^a^/(7r e'^). There is one special case, where there is no need to use any additional assumption like 
Eq. (8), the case of negligible inter-bilayer conductivity. Then Aamt can be neglected and we obtain 



dint 

dl 



A{Hc) — k ^ Aa or equivalently 



"■bl 



AH,[meV] = 4.7 • 10 ' ^fa' + ^^[^-1,^-2] (^0) 
dbi [A] du 

The following points deserve some comments. 

(i) Applicability of Eq. (10): Equation (10) can certainly be used for systems like Bi-2212, where the inter-bilayer 
conductivity is known to be very small. However, it can still be used — with a precision of about 20% — for the less 
anisotropic Y-123 because even there the ratio ^^fnt/'^u obtained from the infrared data is fairly small, especially 
in the underdoped samples. 

(ii) Comparison with the sum-rule for single layer materials: Note the difference between the right hand sides of 
Eq. (3) and Eq. (10). The right hand side of Eq. (10) is larger by a factor of {du + dintY /d\i than the one of Eq. (3) 
(with d = du+dint)- It means that the sum- rule in Eq. (3) underestimates A{Hc) by a factor of df^ / {du + dint)"^ (V^O 
for Bi-2212). For a strongly anisotropic bilayer compound even a tiny change of SW may correspond to a significant 
change of (He). 

(iii) Interband polarizability: In the derivation of Eq. (10) we have tacitly assumed that the interband polarizabilities 
of the two regions, intra-bilayer and inter-bilayer, are the same, i.e., that the two conductivities, au and aint, contain 
the same term —iujEao- For e^oU 7^ ^oo int we would obtain 

A(jj^) =k^=k ^due^^nt+d^ntsc^uf 
du dl{du+d,nt)el,,nt 

(iv) SR for multilayer cuprates with n > 2: The approach presented above can easily be generalized and we obtain 

A{H,) = k (^-1H'+^'"* A^. (12) 
du 

The distance between the closely-spaced copper-oxygen planes is denoted by du, dint is the distance between the 
multilayer blocks. In order to obtain Eq. (12) we have assumed that the inter-multilayer conductivity is negligible and 
that all the "Josephson junctions" (regions between neighboring copper-oxygen planes) within the multilayer block 
exhibit the same electronic conductivity. 

In Table I we show the values of ps, AN, Aa ^ ps + AN, A{Hc), Ej, and the condensation energy per unit cell 
{Uq) for Y123, Bi-2212, and Bi-2223. 

Ps The values of ps for Y-123 have been obtained by using Eq. (A3) of Ref. ||] and the values of uju and ujint presented 
in Refs. [^J|]. The corresponding values of the c-axis plasma frequency are 250, 520, and 1500 cm^^. The values of ps 
in the Bi-compounds are negligibly small. 

AiV AiV = Ns{n,) - Nn{^c), Ns{uj) = N{T « T^cj), 7V„ = N{T w T„uj), N{T,u;) = aic(r, t^') dc^'- The 
values of AA^ for underdoped Y-123 have been obtained by integrating the conductivity data presented in Refs. [^|jl^] 
supplemented by mid-infrared data ranging up to flc = 1500 cm^^ that have been more recently obtained by ellipso- 
metric measurements. A linear extrapolation of ctic below 100 cm^^ has been used. Figure 2(a) shows a complete set 
of ellipsometric data for the Tc = 80 K sample including the extrapolation. Note that ctic(T = 20K << Tc) is larger 
than aic{T = 100 K « Tc) in two different regions: in the region around 550 cm~^ (label B in Fig. 2(a)) and in the 
one located above the frequency range of the apical oxygen modes and centered at 900 cm~^ (label C). The additional 
absorption band B has already been attributed ||,^ to the transverse plasma excitation (TPE) |@-^,|ll| . We propose 
that the band C also belongs to the TPE since a splitting of the SW of the TPE into two parts — one below and one 
above the frequency range of the apical-oxygen modes — is consistent with the Josephson superlattice model (JSM) 
[f7[-p|,pl]|. This is demonstrated in the inset of Fig. 2(a) which shows results of the model calculations of aic for two 
sets of parameters, one corresponding to the SC state and one to the normal state. Details of the calculations are 
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given in Ref. ||T^. The frequency dependence of the quantity Nn{uj) — Ns{uj) corresponding to the data of Fig. 2(a) 
is shown in Fig. 2(b). Note that our value of |A7V| of 0.5ps is considerably smaller than the one of O.Sps obtained by 
Basov et al. ||] for samples with a similar value of T^. This is mainly due to the fact that the band C appears above 
the value of ilc of Ref. (800 cm^^) while below our value of 1500 cm"'^. The value of AiV for optimum doped Y-123 
has been estimated from the inset of Fig. 1 of Ref. |^ . The difference with respect to the result presented in Ref. |^ 
is again related to the difference in the value of the cutoff frequency Vic- The values of AA^ for Bi-2212 and Bi-2223 
are taken from Refs. Q and respectively. 

A(iJc} The values of A(i7c) resuh from Eq. (9) (Y-123), Eq. (10) (Bi-2212), and Eq. (12) (Bi-2223). 

Uq The values of Uq have been obtained from the specific heat data in Ref. |l^ (Y123) and in Ref. ||l^ (Bi-2212), 

the value for Bi-2223 is an estimated upper hmit, ;7o(Bi - 2223) « 2C/o(Bi - 2212). 

In all the compounds studied {He) decreases below Tc and in all of them except for Bi-2212 A{Hc) is of the same 
order of magnitude as Uq. The values of A{Hc) in Y-123 are smaller (ca by a factor of 2) than those of Ej presented 
in Ref. ||^, whereas according to the simplest version of the JSM, neglecting the single particle tunnelling, they should 
be larger by a factor of 4 MM- This is related to the fact that the values of the SW change Aa are considerably 
smaller (ca by a factor of 8) jM than the estimates of the SW of the TPE based on the JSM. We are aware of the 
following reasons for this discrepancy. 

(i) Non zero single particle contribution, i.e., ANint < and ANu < (see Fig. 1 for a definition of ANint and 
AiVf,;). Physically this means that parts of Psint and psu originate from the low- frequency regions of Uinti and f7f,;i, 
respectively. As a consequence, the SR-based estimate of A{Hc) determined by the changes of the total FIR SWs, 
Aaint < Psint and Aati < PsU, (see Eqs. (6a) and (6b)) is smaller than the JCE estimate determined solely by Psint 
and ps bi . This is probably the main reason for the discrepancy. 

(ii) Finite compressibility (FC) effects. The latter have been shown ||ll| to shift both the frequency of the (longitudinal) 
intra-bilayer plasmon and the frequency of the TPE towards higher values with respect to the "bar values" of the 
simplest version of the JSM. At the same time the FC effects do not influence the SW of the TPE. Since the effects 
have not been considered when analyzing the data in Refs. [p|-p^P], the resulting values of ujm presented therein, and 
collected in Ref. H], may be a bit higher than the actual ones. This would imply that also the values of Ej in Ref. j|] 
are higher than the actual ones, which would account at least for a part of the discrepancy. We do not think, however, 
that this is the main reason for it. If the influence of the FC effects was considerable, the frequency of the TPE in 
Bi-2223 would be significantly lower than in Bi-2212 because the FC induced shift of the plasma frequencies decreases 
with increasing distance between the outer copper-oxygen planes of the multilayer block. This has not been observed 

mm- 

(iii) "Pseudogap below Tc" . The difference between the results of the two approaches (JCE and SR) could be much 
smaller, if we considered the temperature evolution of the normal state spectra below Tc, i.e., if Aa was taken as 
a(T << Tc, r2c) — limT_^oCK (normal state, T, flc). It is this limit that should be actually used when estimating A{Hc) 
as emphasized also elsewhere |2]. 

(iv) Fluctuation effects above Tc. In strongly underdoped Y-123 the additional absorption peak starts to form already 
at temperatures much higher than Tc |20 ij], presumably due to fluctuation effects ||,^. The JCE-based estimate of 
A{Hc) contains a contribution of these effects since it is determined by the low temperature value of cobi (i.e., in a 
way that does not require any assumptions concerning the onset of superconductivity) . On the other hand the effects 
obviously do not contribute to the SR-based estimate determined by the change of the spectra below Tc- 

The SW of the additional absorption band corresponding to the TPE (Aa in Table 1) is suprisingly small in 
Bi-2212 as compared to Y-123. We suggest that this is largely due to the fact that the spacing layer in Bi-2212 is 
more insulating than that of Y-123. Indeed, it can be seen from Eq. (11) that for a given value of Aa;,; a decrease 
of the ratio Eooint/ ^ooU leads to a reduction of Aa. A similar reduction of Aa can also be caused by the electronic 
background of the intra-bilayer region: this is what has been assumed in Refs. [p^p|]. These observations indicate that 
the estimate of A{Hc) in Bi-2212 based on Eq. (10) may be considerably lower than the actual value of this quantity. 
Furthermore, it can not be excluded yet that a part of Aa appears at frequencies above the FIR range. At present 
we do not know why Aa is almost by a factor of 5 larger in Bi-2223 than in Bi-2212 (the JSM yields a factor of ca 2). 

In summary, we have developed an extension of the c-axis tight-binding SR that applies to multilayer HTCS and 
allows one to estimate — model independently — the kinetic energy change A{Hc) associated with the SC transition. 
For multilayer HTCS with insulating (or almost insulating) spacing layers the ratio between the SW change and A{Hc) 
is determined by a geometrical factor that is typically an order of magnitude lower than the one of the conventional 
tight-binding SR. Using published far-infrared data that are in part complemented by new MIR data for underdoped 
Y-123, we found that {He) decreases below Tc in all the compounds studied including optimally doped Y-123 and 
almost optimally doped Bi-2212 and Bi-2223. The decrease seems thus not to be restricted to the underdoped regime. 
In all the compounds studied except for Bi-2212 A{Hc) as determined by the SR is of the same order of magnitude 
as the condensation energy and there are several reasons to believe that for Bi-2212 A{Hc) is underestimated. To 
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conclude, in all the multilayer HTCS studied the decrease of {He) below Tc does represent a significant contribution 
to the condensation energy. Its high value suggests that it may well be the dominant contribution. The possibility, 
however, that the changes of {He) below Tc simply parallel much larger changes of the in-plane kinetic energy (such 
as obtained in Ref. |2l) for Bi-2212) cannot be excluded. 

We acknowledge stimulating discussions with D. van der Marel. T. H. thanks the AvH Foundation for support. 
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TABLE I. Values of the superfluid spectral weight ps, the quantity AA'' defined in the text, the spectral weight change Aa, 
the kinetic energy change A{Hc), the Josephson coupling energy Ej as obtained in Ref. |Q, and the condensation energy Uo- 
The quantities ps, AN, and Aa are given in mf2~^cm~^, the quantities A{Hc), Ej, and Uo in meV. 



compound 


T4K] 


ps 


AiV 


-AN/ps 


Aa 


A{H,) 


Ej 


Uo 


A{Hc)/Uo 


YBa2Cu306.55 


53 


1.7 








1.7 


0.08 


0.13 


0.05 


1.6 


YBaaCugOs.Ts 


80 


7.2 


-3.5 


0.5 


3.7 


0.16 


0.30 


0.16 


1.0 


YBaaCuaOe.gs 


91 


60.0 


-48.0 


0.8 


12.0 


0.47 


1.14 


0.36 


1.3 


BiaSraCaCuaOg 


91 





0.3 




0.3 


0.02 


0.13 


0.13 


0.15 


Bi2Sr2Ca2Cu30io 


107 





1.4 




1.4 


0.11 




0.26 


0.4 



FIG. 1. Schematic representation of (a) the bilayer geometry and (b) the spectra of ai„ti = Reio'int} and atn = Re{aM} 
and the quantities that describe the related spectral weight changes (as discussed in the text). 

FIG. 2. (a) Experimental spectra of the c-axis conductivity of slightly underdoped (Tc = 80 K) Y-123. The labels A, B, 
and C indicate the region of a pronounced gap- like suppression of ctic, the main part of the additional absorption band due 
to the transverse plasma excitation, and its high frequency satellite, respectively. Inset: results of model calculations for the 
superconducting state (solid line) and for the normal state (dashed line) demonstrating the splitting of the additional absorption 
band, (b) Spectra of the quantity Nn — Ns defined in the text. 
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Fig. 2 
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